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SUMMARY 


The s imi larity rules for linearized compressible flow theory 
( Gothert ’ s rule and its supersonic counterpart) are extended to second 
order. It is shown that any second-order subsonic flow can be related 
to "nearly incompressible" flow past the same body, which can be 
calculated by the Janzen-Rayleigh method. 


INTRODUCTION 


The linearized small-disturbance theory of steady compressible flow, 
based on the Frandtl-Glauert equation, yields a first approximation for 
thin objects moving at either subsonic or supersonic speeds. More pre- 
cisely, it provides the first term in an asymptotic expansion of the 
solution for small disturbances, provided that the flight Mach number 
is not too close either to unity (transonic flow) or to infinity 
(hypersonic flow) . 

The similarity rule that governs linearized subsonic flow past 
general three-dimensional objects was first given correctly by Gothert 
(ref. l) . It has an obvious counterpart in supersonic flow, and the 
rules have rendered great service in both theoretical and experimental 
investigations . 

Recently, various investigators have sought to improve on the line- 
arized theory by finding higher approximations (see, e.g., refs. 2 to 5 ) . 
The second step is commonly referred to as the second-order small- 
disturbance theory, or simply "second-order theory." It can be found 
in general by iterating upon the linearized solution, retaining all terms 
out of the full nonlinear equations of motion whose contribution is of 
the order of the square of the disturbances in linearized theory (ref. 3 ) * 
In the simplest case of plane flow without stagnation points, the line- 
arized disturbances are proportional to the thickness ratio t of the 
airfoil, so that second-order theory adds terms in t 2 , and higher 
approximations extend the series in powers of t 41 . Stagnation points 
lead to the appearance of logarithmic terms, beginning with r^Zn t in 
the fourth approximation. The series diverges in the immediate vicinity 




2 


NACA TN 3875 


of stagnation points, although it can he corrected there by simple tech- 
niques (ref. 6) . Slender pointed objects, such as a smooth body of 
revolution, cause smaller disturbances than airfoils, but logarithmic 
terms always arise at the outset; hence tfye linearized solution contains 
terms in T 2 7n t and r 2 , and the second- of der increment then consists of 
terms in T 4 ln 2 r, t, and t 4 . Nothing is known of the convergence 

of these series; they are perhaps only asymptotic expansions for small 
thickness. Second-order theory, like linearized theory, breaks down in 
the transonic and hypersonic ranges, though it may penetrate somewhat 
farther into their fringes. 

A similarity rule for second-order theory has recently been given in 
the special case of supersonic flow past ^hin flat wings by Fenain and 
Germain, who demonstrate its usefulness in theoretical studies (ref. 5) • 
However, as in linearized theory, the rules for flat wings are only 
special cases of those for general three-dimensional shapes. The present 
paper is devoted to deducing the general pules for subsonic and supersonic 
flows, and examining their implications. In particular, it iB Bhown how 
the rule for subsonic flow relates the second-order solution for any 
object to nearly incompressible flow past the same body, which can be 
calculated by the Janzen-Rayleigh method. 


DERIVATION OF RULES FOR BODIES OF REVOLUTION 


A body °f revolution is the simplest: shape that is not a special 
case, but displays the full generality of the existing similarity rules 
for subsonic, supersonic, transonic, and hypersonic flows. The same can 
be shown to be true of the second-order rjules to be discussed here. 

Hence for clarity of exposition, the secohd-order rules will be derived 
in detail only for an axi symmetric body at zero angle of attack. The 
rules for general three-dimensional thin or slender objects will there- 
after be stated without proof. The subsonic and supersonic cases are 
so s imi lar that they can be treated simultaneously i 

• - tR(x) , where t is a thickness 
say, the maximum slope, average slope, 
thickness ratio, or the like), and 
r(x) is a function of order unity 
(sketch (a)). As usual in similarity 
analysis, the characteristic slope t 
is regarded as a parameter, bo that 
different values of t correspond to 
affihely related members of the same 
family of bodies. 

To second order the flow is 
irrotational, so that there exists 
a velocity potential $(x,r;M,7,r) . 


Let the body be described by 1 
parameter or characteristic slope (e 



Sketch (a).- Notation for 
body of revolution. 
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This notation indicates that for each family of bodies (associated with 
a given function R(x) ) , the flow field is regarded aB depending not only 
upon the two independent variables x and r but also upon the three 
parameters following the semicolon: 

M free-stream Mach number 

7 adiabatic exponent of gas 1 

r characteristic slope of body 

The aim of a s imi larity analysis is to transform the problem so as 
to reduce the number of parameters appearing in it. If that can be 
accomplished, flows having different values of the original parameters 
are related provided only that the reduced parameters are equal. The 
transformation to be used here consists in separating the dependent vari- 
able $ into several components, and then stretching each component and 
the independent variables by factors that depend upon the original param- 
eters. It Is convenient, and involves no loss of generality, to leave 
streamwise coordinates unchanged, so that r is to be stretched but 
not x. 

Perturbation potentials are first introduced by setting 


-£- = x + 0+ cp + ... (l) 

Uqo 

where 0. is the potential of linearized theory, and cp the second-order 
increment . 

Rules for Linearized Theory 


The linearized problem is 


0z 


□ i ■ (1 - + 0rr + ^ = 0 


0 


at infinity 


( 2 ) 


f T = tR’ (x) at r = tR(x) J 

1 Hayes has pointed out (ref. 7) that to second order an imperfect 
gas corresponds to a polytropic gas having a 7 equal to the free- 
stream value of 


1 



where c is the speed of sound and p the density, the partial 
derivative being taken at constant entropy s. 
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The first relation is the linearized Prahdtl-Glauert equation. The 
second is a statement , sufficiently definite for present purposes, of 
the requirement that the flow approach a ; uniform stream far from the "body 
in alm ost all directions. The third is the linearized condition of tan- 
gent flow at the "body surface. 2 The linearized problem is seen not to 
involve 7, so that the solution depends: upon only the two parameters 
M and t. 

The similarity rules can he obtained by considering arbitrary scale 
transformations of 0 and r. It is readily found that the only choice 
that reduces the number of parameters from two to one is (temporarily 
suppressing the dependence on parameters) 

0 (x,r) = i F(x,p) 

P 

and 

P = 

where 

s/l - M 2 for subsonic flow 
•7 m 2 - 1 for supersonic flow 
Then the problem becomes > 

N 

Fpp + -p- ± ^xx = 0 

F — > 0 at ; infinity 

Fp = £tR 1 (x) at, p = PtR(x) 



(3c) 


(4) 


(3a) 

(3b) 


where here and later the upper and lower signs apply, respectively, to 
the subsonic and supersonic problems. , 

The transformations of 0 and r hhve'been so chosen that the 
problem is reduced to one involving the! tvc_ parameters M and t not 

2 In what is generally called the slender -body approximation, the 
body is assumed to be so smooth and slender that the tangency condition 
can be imposed on the axis rather than bn the actual surface, in this 

case in the form lim r0 r = t 2 R(x)R i (!x) . Thus slender-body theory is 

r->o 

a further approximation within linearized theory (being, in fact, the 
leading term in the asymptotic expansion of the l i nearized solution for 
small thickness t) . Consequently, the slender-body solution obeys the 
similarity rules of linearized theory, and the second-order slender- 
body solution likewise obeys the secon^- order similarity rules. 



NACA TN 3875 


5 


separately, but only in the combination pT. This is the similarity 
parameter. Two subsonic or supersonic flows past bodies of the same 
family are related if the corresponding Mach numbers are such that the 
parameter Pt is the same for both. The nature of the relationship is 
found by reintroducing the dependence on parameters into equation (3 a ) } 
which gives the similarity rules 

0(x,r;M,r) = i F(x,prjp-r) ( 5 ) 

P 


Second-Order Rules 


The second-order problem is found to be (ref. 3) 

□ q> - M 2 |[(r + 1)M 2 + 2(1 - M 2 ) ] 0x0xx + 20r0xr + 0r 2 0rr} 


q> — $► 0 at infinity 

<Pr = T 0x Rt ( x ) at r = tR(x) 


> (6) 


Note that the first equation contains not only quadratic terms on the 
right-hand side, but also the triple product 0 r 2 0 rr "whose contribution 
is of second order in some cases. 

The parameter 7 appears only linearly in the combination (7 + l) 
and can accordingly be separated out . Thus the appropriate transformation 
is found to be 


cp(x,r) 


_ 1 |~f 

' F L 1 


(x,p) + M 2 f 2 (x,p)+(7 + 


1 ) £2 f 3U,p) 




(7) 


Then equating like powers of M 2 yields the following set of three 
problems for f x ,f 2 ,f 3 in which the parameters M, 7 ,t appear again only 
in the form of the single similarity parameter Pt: 


□ fi = 0 

f x — > 0 at infinity > 

f lp = PtF x R* (x) at p = PtR(x) 


(8a) 



6 


NA.CA TN 3875 


□ f 2 - ±2F x Fxx + 2Fp F xp + F p %p 1 

f 2 — > 0 at infinity 

f 2 p =0 at p = |3 tR(x) 

4 

> 

□ f 3 = JxFxx 

f 3 — =► 0 , at infinity > 

f 3 p =0 at p = PtR(x) 


(8b) 


(8c) 


Then reintroducing the explicit dependency on parameters into the func- 
tions f i,f 2 ,f 3 of equation (7) gives thy similarity ruleB for the 
second-order increment in perturbation potential: 

cp(x,rjM,7,T) - -i jf i(x,prjpT) + M 2 f 2 ()+(r + 1) ^ f 3 ()J (9) 

where the arguments of f 2 and f 3 are the same as those of f Hereafter 
the arguments of later terms will be emitted in this fashion when they are 
the same as for the leading term. It should be emphasized that, just as 
in linearized theory, the subsonic and supersonic rules are quite distinct, 
although they have the same form (9). Because of the different defini- 
tions of p, and the resulting ± signs in! equations (4) and (8b), a super- 
sonic flow is not related to a subsonic flow. Discussion of these results 
is deferred to the general case. 


RULES FOR GENERAL. BODIES 




Sketch (b).- Example of two- 
related bodies. 


Consider a family of general 
three'-dimensional bodies, whose mem- 
bers are derived from one another by 
a uniform magnification or reduction 
of all dimensions normal to the free 
stream (sketch (b)). Each member of 
such y family can be characterized, 
as before, by some characteristic 
slope; It may be emphasized that 

t can be identified with thickness, 
camber. Or angle of attack, all of 
which vary together for related 
bodies . 

The preceding analysis can be 
extended in a straightforward way to 
such general bodies, at the expense 
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only of typographic complexity. Both cross -stream dimensions "behave in 
the way that r did before. Hence the subsonic and supersonic second- 
order rules for the velocity potential, corresponding to equations ( 5 ) 
and (9) , are, in Cartesian coordinates: 

jf- ®(x,y,z;M,7,T) - x + i F(x,py,pz;3-r) + \ f i(x,py,pzjpr) + M 2 f 2 () + 

u oo p L 

(r + 1) p UO] (10a) 


Differentiation yields the corresponding rules for velocity 
components (those for w having the same form as for v) : 

- 1 + p' U(x,py,pz; 3 r) + ^ juaXx^y^zjpr) + M®u 2 ()+(7 + l) U 3()J 

(10b) 

^ | V(x,py,pz;pT) + i ^(x^y^z jp-r) + M*v a ()+(7 + D v 3 ()] 

(10c) 

(The functions appealing here are actually related to derivatives of the 
functions in equation (10a) , but the connection is of little interest.) 
To second order the pressure coefficient is given by 


Op - -20x ~ (0y 2 + 0z 2 ) - 2? x “ ^^y^y + 0 z c Pz) + ( m2 “ -^^x 2 + 
M 2 0 x (0y 2 + 0z 2 ) + | M 2 (0y 2 + 0 Z 2 ) 2 


where the terms in the second line may be significant for slender shapes. 
Substituting the expressions (10) for velocity components and simplifying 
shows that the similitude for pressure coefficient has the same form as 
that for the streamwise velocity increment Au/Uo,,: 
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C p (x,y,z;M,r^T) = i P(x,py,pz;pT) + ^ |p.i(x,Py,pz;PT) + M 2 p 2 () + 

fr + D^PaOl (LOd) 

P 

The similarity rules for the perturbation stream function in plane 
flow are the same as those for v/Uoo (eq_. (10c)). 

Alternative Forms. _ 

As with other similarity rules , an unlimited number of alternative 
forms can be produced by multiplying by powers of the similarity param- 
eter. Thus, of the many possible alternatives to the second-order 
rules (lOd) for pressure coefficient, two of the most useful are: 

Cp = | P(x,py,pz;PT) + [p l0 + M^sO+Cy + l) jjg SsOj ( 10e ) 

C p = T 2 P(x,py ,P zjPt) + T^j^PiO + M^sO +(7 + l) ^2 P 3 ()J (lOf) 

In addition, the first two second-order terms can be manipulated, using 
the connection between M 2 and p 2 , to yield additional alternative forms 
such as the following, which correspond to the three forms above: 

G p - i P(x,Py,pz;pT) + ^ Pi() + jjg Vz' ( )+(7 + 1) PsOj ^ 10g ^ 

Cp = ^ P(x,py,pzjpT) + r2^Pi() + ^ i P 2 , ()+(7 + 1) P3()J ( 10h ) 

Cp = T E P(x,Py,pz;PT) + r^p 2 p 1 ( ) + p 2 *()+(7 + l) P3OJ ( 10i ) 



NACA TM 3875 


9 


Force Coefficients 


The rules for pressure imply rules for the lift and drag coefficients. 
For example, equation (lOe) leads to 

C L (M, 7 ,t) = | L(pr) + p + M 2 Z 2 (pr )+( 7 + l) ^ * 3 (Pt) 

(10J) 

Cp(M, 7 ^T) = ^ D(pr) + fdi(PT) + M 2 d 2 (pT)+( 7 + l) ^ ^sCPt) 
p p L p J 

( 10 k) 

if the coefficients are referred to some plan-form area. If some cross- 
sectional area is used, each term is reduced by one power of t. Various 
alternative forms are again useful. In the case of lift coefficient, 
one will ordinarily choose to identify r with the angle of attack. 


RULES FOR QUASI- CYLINDRICAL BODIES 


A special class of objects must be distinguished, which will be 
termed quasi -cylindrical bodies. These are shapes that lie everywhere so 
close to seme cylinder (not necessarily circular) parallel to the free 
stream that to a first approximation the condition of tangent flow can be 
imposed at the cylinder rather than on the actual body surface. Likewise, 
in second-order theory the tangency condition can be transferred to the 
cylinder by Taylor series expansion. The simplest example is an airfoil 
whose thickness, camber, and angle of attack are so small that the tan- 
gency condition can be transferred from the airfoil surface to a mean 
plane parallel to the stream (sketch (c)). Another example is an open- 
nosed body of revolution whose radius varies only slightly. Others are 
biplanes, cruciform wings, any of 
these in an open or closed wind 
tunnel, in combination with one 
another, etc. 



A quasi- cylindrical body can be 
regarded as consisting of a skeleton 
upon which Is superimposed a small 
slope distribution. The skeleton is 
simply the projection of the body 
onto the basic cylinder. For 
example, the skeleton of the quasi- 
cylindrical body of revolution is 
the circular tube Bhown dashed in Sketch (c).- Examples of quasi- 
sketch (c). cylindrical bodies. 



77777777 
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The special place of quasi- cylindrical "bodies in similarity theory 
arises from the fact that the skeleton and the slope distribution can be 

_ — - varied independently. This extra 

fY~ ~ freedom is important. For example, 

|| J j it leads to a useful transonic siml- 

\7 | | ~ larity rule for quasi-cylindrical 

' — — — ~ bodies whereas none exists for general 

shapes. -It is convenient always to 
■ 1 1 leave; streamvri.se dimensions unaltered. 

' \ u Hence;, 'vTe - consider families of quasi- 

| cylindrical bodies that are derived 

I Ij from one another by a lateral com- 

J Ij pressfon or expansion of the skele- 

V/ || fm _ y ton, and_.a quite independent magnifi- 

^ " cation or reduction of all surface 

Sketch (d).- Example of two slopes* Two members of such a family 

related quasi-cylindrical are shown in sketch (d). 
bodies . 


Distortions of the skeleton will be measured by some characteristic 
"aspect ratio" A. It is important to note that the term "aspect ratio” 
is used here in a very general sense to mean any typical ratio of gross 
cross stream to streamwise dimensions. Fbr example, in the last shape 
in sketch (c), the ratio of wind-tunnel height to airfoil chord is an 
appropriate characteristic aspect ratio. * Changes of slope are measured, 
as before, by some, characteristic slope T. 

The preceding similarity rules can be simplified for quasi-cylindrical 
bodies by using the facts that first -order perturbation quantities are 
directly proportional 3 to t, and second-order terms to t 2 , The simpli- 
fication can be carried out by first imagining the quasi-cylindrical body 
to be restricted to be a general body, which means that both pA and |3t 
must be the same for similarity. Then consider the preceding rules in the 
particular* alternative forms In which r \ appears explicitly outside the 
first-order term and t 2 outside the second-order terms. For the 
pressure coefficient, this form is that of equation (lOe): 


C P = f 



+ M 2 p 2 ()+(7 + l) 




So far the functions P,pi,p 2 ,P 3 have been supposed to depend para- 
metrically upon both (SA and |3t. However, the first- and second-order 
terms can be proportional to t and t 2 , respectively, only if the supposed 
dependence upon |3t Is nonexistent. Hen'ce, 1 the similarity parameter 
is pA alone, and the rules for pressure, are. (dropping bars from the 
functional symbols) : 

3 This is by no means true for general bodies; as noted previously, 
the first-order pressure coefficient on a; smooth slender pointed body of 
revolution varies as r 2 Zn t for small t. 
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Cp(x,y,z;M,7,*r,A) = P(x,£y,| 3 zj 0 A) + Px(x, 0 y, 0 z;pA) + M^aO + 

p p L 


(7 + 1) 


p - 0 ] 


(lla) 


The corresponding rules for the potential and velocity components can, 
if desired, be written down by inspection from equations (10) . 

With p 2 = P2 - Pi (where, as before, the upper sign applies to 
subsonic and the lower to supersonic flow) , these rules can be rewritten 
as 

Cp = ^ ?() + t2 |p:lO + ~a ^2(^ + ^ + l) ^4 Pa() J (lib) 

and this is the result that Fenain and Germain found in their treatment 
of the flat diamond cone in supersonic flow (ref. 5). 


Connection With Hayes' Rule 


For plane flow past a single body, Hayes has discovered a remarkable 
rule for the second-order surface pressure (ref. 7) • It implies that, 
on the surface, the functions in equation (lla) are such that p 2 = 0 
and Px = %> 3 . Hence, 


c Pe . 1 PM 


Ki - m 2 )- 


Pi(x) 


(12) 


In supersonic flow this is simply Busemann' s well-known second-order 
solution, P being twice the local slope of the surface and p x twice 
its square; in subsonic flow P and px are more complicated (ref. 8) . 

This rule implies a corresponding, but more complicated, rule for 
surface velocity (ref. 8). 

In addition to the restriction to single bodies and plane flow, 
these rules are not similarity rules in the sense of the preceding results, 
because they apply only at the surface rather than throughout the field. 
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EXAMPLES 


The rules will be illustrated by two' simple examples, attention 
being confined to the surface pressure coefficient. 

Slender Circular Cone in Supersonic Flow 

Broderick has derived the second-order slender-body solution for a 
circular cone at zero angle in a supersonic stream (ref. 2 ) . The surface 
pressure coefficient on a cone of slope f is 

Cpg = ^(2 in - 1 ) + T«|W*(ln A) - (5M? - l)2n ^ M M 2 + i + 

(r + 1) (13) 

This has the form of equation (lOi) with ' 

S = 2 m If - 1 : p»* - 111 h * ¥ 

^ - 3 in * r, - ’ # - f . > - 1 

Wavy Wall in Closed Subsohic Wind Tunnel 


Consider the sinusoidal wall 
y = t sin x at a distance h from 
a flat wall (or a distance 2 h from 
its mirror- image) as indicated in 
sketch (e). Subsonic flow between 
the walls at a mean Mach number M 
can be readily calculated to second 
order by separation. of variables. 

The resulting pressure coefficient 
on the surface of the wavy wall is 


N /'~S X N 

\ / \ / \ 
s_ ^ s.y’ 

y 

1 hj/ / // // // / / 111 


1 



Sketch (e) .- Wavy wall in 
wind tunnel. 
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2rl±e^ (7 + l)M 4 +lt(l-M 2 ) 1 + 4e~ 2f3h + e“ 41311 

* 1-e- 2 ^ 4(l - M 2 ) 2 (IT^P 


2T 2 ^ +2(7 + 1) ^ T 2 P— 5 ~y 5 COB 2x (l 4 ) 

(l-e“^) P (l + e- 2 ^)(l-e-^) 


The relevant aspect ratio is the height h (which is really a multiple 
of the height-chord ratio, because of the choice of scale for the wavy 
wall) . Thus the result is seen to have the similitude of equation (lla) . 

As the tunnel height increases indefinitely, the last two terms disappear, 
and the remainder follows the similitude of equation (12) for the surface 
of a single plane body. 


REDUCTION OF SUBSONIC PROBLEM TO 
NEARLY INCOMPRESSIBLE FLOW 


In linearized theory, an important application of the similitude is 
Gothert’s rule, which reduces any subsonic flow problem to a related 
incompressible flow (ref. l) . As the rule is usually stated, the incom- 
pressible flow is that past a thinner affinely related body. However, 
the incompressible solution for one member of an affinely related family 
of bodies determines that for all other members, so that the subsonic 
flow may, if desired, be related to the same body rather than a thinner 
one, and that viewpoint will be adopted here as being the simplest. 

In second-order theory, the explicit appearance of terms in M 2 
and (7 + l)M 4 in equations (10) means that reduction to an incompressible 
problem is impossible (except for the special case of the surface of a 
single plane body, where eq. (12) applies) . The second-order problem 
can, however, be reduced to a nearly incompressible flow. 

Flows at low Mach numbers can be calculated by the Janzen -Rayleigh 
method, which involves iterating on the incompressible solution to obtain 
a power series in M 2 . Thus the velocity potential is approximated by 

$(x,y,z $M, 7 ,t) = S>oU*yj z ; T ) +M 2 $ 1 (x,y,zj+) + 

u c» 


( 7 +l)M 4 $ 2 ()+ M*tt 3 () + 


( 15 a) 
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The two terms in M 4 are ordinarily considered together, "but for present 
purposes it is essential to separate them because only $> 2 is required. 
This is fortunate because $ 2 can be calculated almost as easily as 
■whereas the determination of $> 3 is much more difficult. 

The small- disturbance and Janzen-Rayleigh series represent two dif- 
ferent asymptotic expansions of the actual solution.. They are believed 
to complement each other, so that an expansion of the Janzen-Rayleigh 
solution for small thickness must be identical with the expansion of the 
small-disturbance solution in powers of M 2 , as has been verified in all ., 
worked examples. This fact permits the sniall-disturbance solution to be 
recovered from the Janzen-Rayleigh series. The converse is not true, 
however, except for bodies without stagnation points, because the small- 
disturbance expansion is not uniformly valid near a stagnation point. 


Procedure for Recovering Second-Order Solution 


Suppose that the Janzen-Rayleigh series of equation (15a) is known 
up to the term in (7 + l)M 4 . Then expanding each term formally for small 
thickness t, and retaining only second-order quantities, gives a double 
series of the form 


— - 3>(x,y,z;M,7,r) = [x + T2$> 01 (x,y,z;T) + t4® 02 () +...]+ 

u oa 

M 2 [ T 2 ® 11 ( ) + 2,2 ( ) + »..] + 


(7 + 1)M 4 [ t 4® 22 ( ) +...]+ o(m 4 ,t s ) (15b) 

Here i|) nl and $ n2 are, respectively, the first- and second-order com- 
ponents of $ n . How also, according to the similarity rule of equa- 
tion ( 10 a) in the alternative form corresponding to equation (lOf), the 
second-order small-disturbance solution hag the form 


77 ~ ®(x,y,z;M,7,T) = x + T^x^y^z j(3t) + t 4 g x () + M 2 gg() + 
u°° l 

(7 + 1) p s 3 ()j + OCt 8 ) (16a) 

and this can be expanded fully in powers of M 2 by using Taylor series 
expansions such as ; 
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s(py) = g(y) - M^g* (y) + o(m 4 ) 

to give 

i~*(x,y, Z ;M,r,T) = x + t 2 -|G(x,y,z;T) - § M 2 jyG y () + zG z () + tG t () J + 

t 4 -| S l() - §M 2 |yg ly () + zg lz () + Tg lT () + 

M 2 g 2 ()+(7 + l)M 4 g 3 ()j» + 0 ( t 6 ,M 4 ) ( 16 b) 

If the two double series expansions (15b) and (l6b) are identical, like 
terms can be equated to give 

G(x,y,z;r) = <Soi(x,y,z j t) 

Si(x,y,z;x) = ( &o2( x ,y, 2 ;T) 

g2(x,y,zjr) = ®i2(x,y,zjT) + ~ Jy$ O 2 y 0 + z$o2 z () + t 0> O2t ()J 

ff3( x ^y^z;T) = $2s(x,y,z;T) 

Then substituting these expressions into equation (l6a) expresses the 
second-order potential in terms of the Janzen -Rayleigh solution as 

jjr $(x,y,z ;M,7 ,t) = x + T 2 $ 01 (x,py,pzjpT) + T 4 -|$ 02 (x,Py,pz;pT) + 

4-0 + ^ *° 2 py () + IT $oa pz () + ¥ ®° 2 pr () + 

+ ® 22 ()j- ( 17 ) 

Here $o2p y ( x >Py>!3 2 SPT) means (<V8y)$ 02 (x,y,z;T) evaluated at x = x, 
y = py> 2 = pz, and t = p T , for example. Note that ® lx is not required. 
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Application to Parabola 


y 



Sketch (f).- Parabola at zero 
angle of attack. 


As an example, consider plane 
subsohic flow at zero angle of attack 
past the parabola described by 

y = tVx (sketch (f)). The Janzen- 
Rayleigh solution including terms 
in (7 + l)M* has been calculated 
by Imai (ref. 9 )- Although the 
velocity potential is complicated, 
it simplifies when expanded in powers 
of t to 


= x + rn + i M 2 *rn 2 ^ ~ i M 2 T 2 ln - — 


U„ 


TTJ 


-TT] 


J; 2 + T) 2 8 


LJlI. m 4 t 2 
32 


i (i 2 + T\ z ) 


+ In 


2\2 r 


i z + T] g 




(18a) 


Here | ,tj are parabolic coordinates related to the Cartesian coordinates 

"by 


[(X - ^ t2 ) + iyj = (I + in)‘ 


so that to second order 


2£ 2 = n/x 2 + y 2 + x I 
2 x\ z = n/x 2 + y 2 - x 


(18b) 


In this case the terms involved in equations (15b) and (17) are 
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IT 


$01 ~ •?: 


$02 — 0 


■ 12 


$22 — 


i_ [1 

8-r 2 L 




I + 


+ In 


juei 

T 2 J 


i^L U 2 + n 2 ) 2 


+ In 




(18c) 


This example illustrates the fact that for planar systems these terms 
are not of order unity in t. Then according to equation ( 17 ) the 
second-order small-disturbance solution is 



x +|n 


1 M 2 
8 p 5 


t 2 ( 4 



+ Zn 



where 


7 + 1 
32 p* 


r 2 k 
— 


i4 


(I 


=2 


+ n 2 )‘ 


+ In 


r . + j 2 

P 2 t 2 




2g 2 = *7x 2 + P^ 2 + x 
2q 2 = Vx 2 + p 2 y 2 - x 


(19a) 


(19b) 


This result is of interest because it apparently cannot be found directly. 
Plane small-disturbance flows can be calculated easily if one adopts the 
thin-airfoll approximation of transferring the boundary conditions to the 
line y = 0 by Taylor series expansion, but that approximation fails near 
round noses in the second approx ima te on and, as a consequence, divergent 
integrals arise (ref. 8). Instead, one can try to treat the round nose 
more carefully using conformal mapping (cf. ref. 11, pp. 361 - 367 )* 
the result is found to be indeterminate to the extent of a multiple of 
Zn(| 2 + q 2 ) . This is the potential of a point source at the origin, 
which is an eigensolution, the proper multiple of which (appearing in 
eq. ( 19 a)) is not determined by the suggested method. 

The second-order increment in equation (l9a) is seen to include 
terms in i^ln t, whose function is to render the argument of the 
logarithm dimensionless. However, these terms are simply constants. 
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so that no logarithms of thickness appeal in the actual flow quantities 
such as velocity and pressure. As remarked~in the introduction, loga- 
rithmic terms in thickness arise in the actual flew disturbances only- 
in the fourth approximation. 

The second-order small-disturbance solution for the stream function 
can in the same way be extracted from ImAi’s Janzen-Rayleigh solution, 
and the result is found to agree with that calculated directly by Kaplan 
(ref. 10) using conformal mapping. It contains no terms in In t. (The 
direct approach succeeds for the stream function, although it fails for 
the velocity potential, because the tang^ncy condition iB imposed on the 
mass flux, which is affected by the above eigensolution. ) Then using the 
connections between the stream function and velocity potential, one can 
verify the correctness of equations ( 19 ) • 


CONCLUDING REMARKS 
Utility of the Rules 


The second-order rules are scarcely suited for correlating experimen- 
tal data, since tests on four related bodieB would be needed in order to 
isolate the four functions involved. That they are, however, useful in 


theoretical analyses has already 





:n pointed out by Fenain and Germain 
in the special case of supersonic flow 
past ; flat wings (ref. 5) • Previous 
investigators had calculated (erro- 
neously, as it turns out) the second- 
order solution for the flat diamond 
cone ! shown in sketch (g) , and carried 
out numerical computations for three 
different Mach numbers and four values 
of the parameter (3 tan A (reported 
in ref. ll) . Because the latter is 
the similarity parameter |3A of equa- 
tionfe (ll), failure to take advantage 
of the similitude resulted in three- 
fold; uniiecessary duplication of 
competing labor. 


The reduction to nearly Incompressible flow assumes importance for 
bodies with stagnation points. The small- disturbance assumption is 
violated, and, as was noted In the example of the parabola, the second- 
order solution consequently cannot be found, directly. For bodies of 
revolution the difficulties appear to be, even more severe.. In such 
cases It Is convenient to calculate the 1 smzen-Ray leigh solution, and 
from it extract the true second-order solution by the procedure outlined 
above . 
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Further Extensions 


The similarity rules can readily be extended to third and higher 
order in the same fashion (except for complications in supersonic flow 
because of the ultimate appearance of significant vorticity engendered 
by curved shod waves) . The similarity parameter remains unchanged; 
the complexity arising in a proliferation of functions multiplied by 

powers of (7 + l^M 211 ^ 2 ^ ) . Likewise, the small-disturbance solution to 
any order can be recovered from the nearly incompressible solution 
provided by an appropriate number of terms of the Janzen-Rayleigh 
solution. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Oct. 18, 1956 


REFERENCES 


1. Gothert, B. H.: Plane and Three-Dimensional Flow at High Subsonic 

Speeds (Extension of the Prandtl Rule). NACA TM 1105, 19 46. 

2. Broderid, J. B.: Supersonic Flow Round Pointed Bodies of Revolution 

Quart. Jour. Mech. and Appl. Math., vol. 2, pt. 1, Mar. 1949, 

pp. 98-120. 

3. Van Dyke, Milton D.: A Study of Second-Order Supersonic Flow Theory. 

NACA Rep. 1081, 1952. 

4-. Fell, J., and Leslie, D. C. M.: Second-Order Methods in Inviscid 

Supersonic Theory. Quart. Jour. Mech. and Appl. Math., vol. 8, 
pt. 3, Sept. 1955, PP. 257-265. 

5. Fenain, Maurice, and Germain, Paul: Sur la resolution de 1’ equation 

regissant, en seconde approximation, les ecoulements d'un fluide 
autour d* obstacles tridimensionnels . Comptes Rendus, vol. 24-1, 
no. 3, 18 July, 1955, PP* 276-278. 

6. Van Dyke, Milton D.: Subsonic Edges in Thin-Wing and Slender-Body 

Theory. NACA TN 334-3, 1954. 

7. Hayes, Wallace D.: Second-Order Pressure Law for Two-Dimensional 

Compressible Flow. Jour. Aero. Sci., vol. 22, no. 4-, Apr. 1955, 
pp. 284-286. 

8. Van Dyke, Milton D.: Second-Order Subsonic Airfoil- Section Theory 

and Its Practical Application. NACA TN 3390, 1955* 



20 


NACA TN 3875 


9. Ttwr -T , Isao: Application of the M 2 -Expansion Method to the Subsonic 

Flow of a Compressible Fluid Past a Parabolic Cylinder. Proc. 

1st Japan Nat. Cong, for Appl. Math., 1952, pp. 3^9~352. 

10. Kaplan, Carl: On the Simll-Disturbarice Iteration Method for the 

Flow of a Compressible Fluid With Application to a Parabolic . 
Cylinder. NACA TN 3318, 1955* 

11. Lighthill, M. J.: Higher Approximations. Section E of General 

Theory of High Speed Aerodynamics, vol . VI of High-Speed 
Aerodynamics and Jet Propulsion, Princeton Univ. Press, 195^> 
pp. 48ii~486. ; 


NACA - Langley Field, Va. 



